The higher-order C„ dispersion coefficients for hydrogen. 
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The complete set of 2nd, 3rd and 4th order van der Waals coefficients, C„ up to n = 32 for the 
H(f s)-H(f s) dimer are computed using pseudo-states to evaluate the appropriate sum rules. A study 
of the convergence pattern for n < f6 indicates that all the Cn<i6 coefficients are accurate to f3 
significant digits. The relative size of the 4th-order to the 2nd-order C^^ coefficients is seen to 
increase as n increases and at n = 32 the 4th-order term is actually larger. 
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The long-range interaction between two spherically 
symmetric atoms can be written in the general form 



V{R -> oo) = -V(,{R) - VsiR) - VioiR) -Vn{R)~..., 

(1) 

where the dispersion potential, Vn(R) of order n is writ- 
ten 



(2) 



The Cn parameters are the van der Waals dispersion co- 
efficients. The even (n = 6,8, . . .) dispersion coefficients 
are calculated using sum rules derived from 2nd-order 
perturbation theory and provide an attractive interac- 
tion. The odd (n — 11,13,...) terms come from 3rd- 
order perturbation theory, and are repulsive IH • Con- 
tributions from 4th-order perturbation theory start at 
n= 12 0,0. 

The dispersion interaction for the simplest system, 
namely the hydrogen dimer, is only known precisely to 
n = 11. The latest calculations by Yan and Dalgarno 
(YD) reported almost exact values for the 2nd-order 
even dispersion parameters up to n = 16. They also gave 
almost exact results for the 3rd-order coefficients, up to 
n = 15. However, it is known that contributions from 
4th-order perturbation theory start at n = 12 Q, H, Q 
so the dispersion potential computed from the YD C„ is 
incomplete for n > 11. 

In this article, the complete set of dispersion param- 
eters, from Cg to C'le are computed by using a large 
pseudo-state expansion to evaluate the appropriate sum- 
rules. The contributions from 4th-order perturbation 
theory to C12, C14 and C'le are explicitly included. 
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All the dispersion coefficients computed in this paper 
were computed by first diagonalizing the semi-empirical 
Hamiltonian in a large basis of Laguerre type orbitals 
(LTOs) defined by 



where the normalization constant is 



i(2Aar) 



(3) 



No. = 



D! 



(4) 



The function L^^^^fLi(2XaT) is an associated Laguerre 
polynomial that can be defined in terms of a confluent 
hypergeometric function Q as 



L^n^tfLi(2Xar) — 



(n„ +£ + 1)! 
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(n„ - l)!(2£ + 2)! 
X Af(-(n„ l),2^ + 2,2A„r) (5) 

All the matrix elements can be written in simple com- 
pact forms provided a common Aq is adopted throughout 
the calculation. However, in the present work, the ra- 
dial wave functions were placed on a numerical grid and 
all matrix elements computed by gaussian quadratures. 
This was done for reasons on convenience as the diago- 
nalisation could be done with an existing program used 
in previous calculations of the dispersion parameters and 
the structures of positronic atoms • This program 
can achieve close to machine precision in almost all radial 
matrix computations. 

Once the Hamiltonian diagonalisation is complete, sum 
rules involving radial matrix elements were used to deter- 
mine the dispersion parameters. The specific sum-rules 
used are those derived by Ovsiannikov, Guilyarovski and 
Lopatko (OGL) 1]. Their expressions are a bit simpler 
than those developed by other authors 0, . There were 
some omissions in their published equations Q, and a 
more thorough description of the sum-rules is presented 
here. 
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TABLE I: The dispersion coefficients for the H-H dimer. All the results in the "Best Previous" column come from [3l except 
that for Cjj' which is taken from Q. All values are in atomic units. 



Coefficient 


10 LTOs 


15 LTOs 


20 LTOs 


Best previous 


Ce 


6.499 026 705 3305 


6.499 026 705 4057 


6.499 026 705 4059 


6.499 026 705 4058 


Cs 


124.399 083 58236 


124.399 083 58362 


124.399 083 58362 


124.399 083 58362 


Cio 


3285.828 414 9425 


3285.828 414 9674 


3285.828 414 9674 


3285.828 414 9674 




1.214 860 208 9619x10^ 


1.214 860 208 9686x10^ 


1.214 860 208 9686x10^ 


1.214 860 208 9686x10^5 


r<-2) 


6.060 772 689 1671x10'^ 


6.060 772 689 1917x10^* 


6.060 772 689 1917x10'' 


6.060 772 689 1921x10® 




3.937 506 393 9865x10** 


3.937 506 393 9985x10** 


3.937 506 393 9985x10* 


3.937 506 393 9992x10** 


Cn 


-3474.898 037 8919 


-3474.898 037 8822 


-3474.898 037 8822 


-3474.898 037 8822 


Cl3 


-3.269 869 240 4549x10^ 


-3.26 986 924 04407x10^ 


-3.26 986 924 04407x10^ 


-3.26 986 924 04407x10^^ 


Cl5 


-2.839 558 063 3179x10^ 


-2.839 558 063 2998x10^ 


-2.839 558 063 2997x10^ 


-2.839 558 063 2998x10'' 




1241.587 803 8317 


1241.587 803 8462 


1241.587 803 8462 


1241.588 




3.009 633 558 9570 xlO^ 


3.009 633 559 0035 xlO^ 


3.096 633 559 0035 xlO^ 




'-'16 


4.745 455 287 4168x10^ 


4.745 455 287 4083x10^ 


4.745 455 287 4079 xlO'^ 




Cl2 


1.227 276 087 0002x10^ 


1.227 276 087 0071x10^ 


1.227 276 087 0071x10^ 


1.227 27609 xlO'^ " 


Ci4 


6.361 736 045 0628x10*^ 


6.361 736 045 0920x10^* 


6.361 736 045 0921 xlO** 




Cl6 


4.412 051 922 7282x10** 


4.412 051 922 7393x10** 


4.412 051 922 7393x10* 





"This entry adds to BM c[^^ to the YD cjj^. 



A. The 2nd-order terms 

The 2nd-order dispersion coefficients for the H-H sys- 
tem have been determined to high accuracy even 
for high n. The working expression adopted for compu- 
tation is, 



^(2) _ 



A+1 

E 



(2A + 4)! 



E 



(2^1 -f l)!(2n -f 1)! 
(E,, + Ei,^ - 2Eo) 



(6) 



where i[ = A + 2 — £i. The state vector repre- 
sents the radial part of the state ii with orbital angular 
momentum ii and energy Ei^ . The ground state energy 
is Eq. The sum rule 



TW=J2{0,0\r%ir 



{2e + 2)i 

2(2^+1) 



(7) 



is a useful diagnostic check of the accuracy of the under- 
lying discetization of the H-spectrum. 



B. The 3rd-order terms; Cn and C13 



The dispersion coefficients, CVi and Cj3, arise from 
3rd-order perturbation theory 0, □,a,y,UJ Close to ex- 
act dispersion parameters for the H-H system have been 
published 



The general expression for the 3rd-order C2A+11 is Q 

C2A+11 = - E E E 

£lki£2 I'-^k'-^l'^ 2li^22^2 

G(A,^i,fi,4,4,A:i,A:i) 



X (o,o|/Mn,^i>(n,^ik''M*2,^2) 

X (i2,^2|r^lO,0)(0,0|/'iK'i,fi) 
X (zi,^'i|r'='H*^,4>(*2,4k'H0,0) 



(8) 



with the notation i?. 



E, 



vjj^i-^ — xjjj I Ei'^ being used in the 
energy denominator. The parameter A is defined 



2A + 8 = ^1 + /ci + £2 + + fci + ^2 



(9) 



and all of the angular momentum indices are greater than 
zero. Defining J = {e^ + ki + 12)12 and J' = {l\ + /sj + 
t'-2)l'2, the coefficient G is defined as 



G(A,4,^'i,^2,4,^i,fci) = (A-H4)!A(J,£i,fci,£2) 
X A( J', £'1, fci , 4)B(A, h^\)B{\ h,k[)Bi\, £2,4) 



where 

and 

A{J,h,ki,£2 



[2(A + 4-4-^2)]! 
(A + 4-£i - ^2)! 



J! 



(10) 



(11) 



(2J + l)!(J-£i)!(J-fci)!(J-^2)!' 

(12) 



3 



C. The 4th-order contributions to C12 and C14 

As far as we know, there have only been two exphcit 
calculations of the 4th-order contribution to Ci2- Bukta 
and Meath fsj gave estimates of cj^ and c['^ for the hy- 
drogen dimer. Ovsiannikov et al [ij developed a general 

and compact expression for the evaluation of Cn^"* , and 
in addition they reported values of 0^2' for all possible 
combinations of hydrogen and the alkali atoms. Rectify- 
ing some omissions in their published equations one 
writes 



'-'2A+12 - 02A+12 - l^2\i 



+6 



a2A2+6 : 



Ai=0 



where 



A = Al + A2 



(13) 



(14) 



The factor a2A2+6 is 



a2A2+6 



A2+1 

E 



(2A2+4)! 



Z_.^ (24 + l)K2£[ + 1)! 



where 



(15) 



(16) 



The expression for a2A+6 is practically the same as eq. 
(2) 

for C^x^f^ ; the only difference being an extra factor in the 
energy denominator (compare with eq. (10) of pj). 
The factor 62A+12 is more complicated and defined as 

Wr2= E E E E E 

^li^2:^3 l'i.l'2:^':^ ^1,^2,^1:^2 ^ 212^22*2*3*3 



Li\Ki\K2\L:i\ 



2£i!2fi!2fci!2fc;!2fc2!2fc^!24!2^!5! 

X (40fciO|£20)(fiOfc;OK20) 
X (fc2040|40}(fc2040|40) 
X {LiQKiQ\m){K2QL30\KQ) 



1/2 



X < 



fci k[ Ki 

f 



h K 



£2 I'o K 



(Ei^i'^ - 2Eo){Ei2i'^ — 2Eo){Ei^i'^ 
X (0,0|/i|n,£i}(ii,4|r'=M*2,^2) 

X (l2,£2|r'=^N3,4)(«3,^3|r'lO,0) 

X (0,0|/M*'i,^'i)(*'i,4|r'=i|*'2,4) 
X (*'2,^2k'n4,4)(4,^3k'H0,0) 
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where Li = (i + £[, L3 = £3 + C^, Ki = ki + k[ and 
K2 = k2 + k'2. We use L = (2i + 1). The sums are 
constrained by the condition 



Li+Ki + K2+L3 = 2X 



(18) 



While £1, £[, £3, £3 must be greater than 0, it is possible 
for £2 and £'2 to be equal to 0. None of fci, k[, ^2 or ^2 
can be zero. Since £2 and £'2 can both be equal to zero, 
the possibility of 12,12 both occupying the ground state 
must be explicitly excluded from the summation. 



III. RESULTS OF THE CALCULATION 

The results of the calculations for the complete set of 
dispersion coefficients up to Cie are given in Table 
The parameters are given for basis sets with 10, 15 and 
20 basis functions per angular momentum respectively. 
The exponent in the LTO was chosen to be A = 1.0 for 
all angular momenta. This choice resulted in much faster 
convergence of the dispersion parameters than that ob- 
served by Yan and Dalgarno in their calculations of the 
3rd-order dispersion coefficients. Table ^ also gives re- 
sults reported by YD and a single calculation of (7^2^ by 



Bukta and Meath (BM) 3 



The 2nd-order terms 



(2) 

The calculations of Ck do not give new information 
and Yan and Dalgarno have given values which are 
converged to better than 15 significant figures. The 
present calculations with the iV = 20 basis are identical 
to 13 significant figures. The small differences in the last 
digit for some of coefficients arise from minor inaccura- 
cies with the radial matrix elements. Hence we conclude 
that the present calculations are numerically reliable and 
that the pseudo-state representation of the H-spectrum 
is close to converged. 

Besides the dispersion coefficients, the sum-rule, eq. Q 
was evaluated and seen to be correct to 12 significant 
digits for all polarities relevant to the evaluation of Cq to 

•-"16 • 



(17) 



B. The 3rd-order terms 

(3) 

Since the 3rd-order terms, Cn have already been given 
by YD, these calculations merely serve as a test of our 
numerical procedures. Once again, calculations with the 
20 LTO basis agree with the YD results to 14 significant 
figures. It is worth noting the present results required 
fewer terms than YD to achieve convergence. YD made 
the choice Aq = 1/(^+1) in eq. © and did not achieve 
convergence to the 14th digit place until the dimension 
of the LTO expansion was 50. The present basis with 
Aq = 1.0 achieves the same level of convergence with 20 
LTOs. 
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TABLE 11: The n > 16 dispersion coefficients for the dimer. All values are in atomic units. 

r^(2) ^(4) ^(2) ^(4) ^(3) 

17 2.726 099 889x10^ 

18 3.234 218 716x10^° 7.009 061 179x10^ 3.935 124 834x10^° 

19 3.020 900 833x10" 

20 3.278 573 440x10^^ 1.083 922 188x10^^ 4.362 495 628x10^^ 

21 3.900 227 980x10" 

22 4.021 082 848x10" 1.832 218 347x10" 5.853 301 195x10" 

23 5.856 636 712x10^^ 

24 5.868 996 335x10^^ 3.444 924 821x10^^ 9.313 921 156x10^^ 

25 1.017 059 252x10" 

26 1.005 294 993x10" 7.249 737 286x10" 1.730 268 722x10" 

27 2.028 440 001x10^° 

28 1.996 944 941 xlO^^ 1.709 243 726x10^^ 3.706 188 667x10^^ 

29 4.613 037 362x10^^ 

30 4.553 288 866 xlO^^ 4.507 006 859x10^^ 9.060 295 725x10^^ 

31 1.188 007 684x10^^ 

32 1.181 107 088x10^® 1.325 398 446x10^® 2.506 505 534x10^® 



C. The 4th-order terms 

The only previous explicit calculation of a 4th-order 
term was that made by Bukta and Meath (BM) Q , and 
the only parameter given was C[t^. TheOGL Q estimate 

of C^j'', 1.220x10^ au was made using an approximation 
to the Greens function and so perfect agreement is not 
expected. However, the present calculation agrees with 
BM calculation ofC['^ to all digits quoted, namely seven. 

The number of terms in the sum, eq. H17II increases 
rapidly as n increases. There are 4 terms for 0^2' i there 
are 64 terms for , and finally there are 460 terms for 

'-^le ■ 

The dominant contribution to C^"^ comes from 62A+12 
with 96% of coming from bi2- The tendency for 

&2A+12 to be the dominant term in C^'^ becomes more 
accentuated as n increases and the bie term gives an es- 
timate to c[q' that is correct to 0.1%. 

One feature of of Table ^ concerns the relative size of 
dn'^ to dn \ For n = 12, the ci^^ : ci^^ ratio is 1.02%. 
However, as n gets larger, the ratio also gets larger. For 
n = 14 the ratio is 4.97%, while for n = 16 the ratio is 
12.1%. 

D. The dispersion coefficients for n > 17 

Higher order contributions than 4th-order begin at 

n = 17. There is a 5th-order contribution to C17 and 

(2) 

a 6th-order contribution to Cis JL]- Estimates of C„ ' 
for n > 17 have been made by a variety of authors 
UEIIlllllIi- However, the only estimate of the 3rd 



and 4th-order terms with n > 17 are those of OGL 0]. By 
explicit calculation they obtained c[^^ = —2.739 x 10^ au 
which agrees with the present more extensive calculation 
to within 1%. Making an approximation to the greens 
function they estimated c['^ = 3.3 x 10^ au, which is 
about half the size of the present value. 

The dispersion parameters up to C32 from the present 
calculation are tabulated in Table ^ The reason for 
taking the calculations so far rests in the relative size of 
the C^^ and Cn^ terms. It was noticed that the C^'' ■.Cn'' 
ratio got larger as n increased. So the calculations were 
extended to C32 in order to demonstrate explicitly that 
the C^^:Cn'^ ratio can actually become larger than 1.0. 

The precision of the entries in Table ^ is not as high 
as those in Table 13 The calculations of clf' did become 
more time consuming as n increased. There were 922,064 
different (£1, /ci, ^2, • ■ ■) combinations by the time n — 
32 was reached. Also the number of radial integrals in 
eq. ()17|l increases as where N is the numbers of LTOs 
for any given I. So the = 20 calculation is 64 times 
more intensive than the A^ = 10 calculation. 

The C4^^ and &n'^ entries in Table HD were taken from 
the calculation with 15 LTOs. The clx^ entries were 
taken from a A^ = 15 calculation up to n = 20, thereafter 
the A^ = 10 basis was used. The values of C„ ' agree with 
those of Thakkar for all ten digits given in Table 
Comparisons between A'^ = 10 and A^ = 15 calculations 
for Cn^ suggest that the convergence is slower as n in- 
creases and that C^^ is reliable to about 6 digits. A 
similar level of accuracy can be expected for and a 
comparison between the A^ = 10 and A^ = 15 values for 
C20'' gives agreement for the first 9 digits. 
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IV. CONCLUSIONS 

The higher n dispersion parameters from Cn, through 
to Ci6 have been computed to an accuracy of 13 sig- 
nificant figures for the H-H dimer. Since the 4th-order 
contributions were included for C12, C14 and Cig the 
adiabatic dispersion interaction can now be regarded as 
complete up to terms of order R^^^. 

The time taken to evaluate the dispersion coefficients 
was not excessive. For example, a calculation using 20 
LTOs took about 17 minutes to determine all terms up 
C16 on a 850 MHz CPU. Hence the pseudo-state method 
adopted here, and in other similar works (e.g. |3i 0|)i 
could be used to make explicit calculations of the 5th- 
order correction to C17 and even the 6th-order correc- 



tion to C18 Q. Therefore, it is certainly possible with 
existing technology to determine the complete dispersion 
interaction for the H-H interaction for all terms up to 
and including C22. 
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